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ABSTRACT 



Modern Naval weapon and sensor systems are strongly 
influenced by the marine environment. Foremost among the 
atmospheric effects is ducting of electromagnetic energy by 
refractive layers in the atmosphere. To assess the effect 
of ducting or. electromagnetic emissions, the Navy developed 
the Integrated Refractive Effects Prediction System (IF EPS). 
Research at Naval Pcstgraduate School (NPS) has led to 
development of a state-of-the-art model which car. be used to 
predict changes to the refractive profile of the lower 
atmosphere. The model uses radiosonde data and surface 
meteorological observations to predict changes ir. refractive 
conditions ar.d lew level cloud/fog formation over. 18 to 30 
hour periods. The model shows some skill in forecasting 
duct regions wh e n subsidence rates can be specified ^o 
within +/-.0015 m/s. This thesis shows the applicability of 
the NPS marine atmospheric mixed layer modal ~o fleet tac- 
tics. Atmospheric refractive effects or; specific emitters 
can be predicted when model predictions are used in conjunc- 
tion with IREPS. 
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INTROD UCTION 



Over the centuries military leaders have learned to 
appreciate and take advantage of the effects cf the environ- 
ment, and when they failed to io so, have met with disaster. 
As it approached the islands of Japan in the thirteenth cen- 
tury, Ghengis Khan's Mongol invasion fleet was destroyed by 
a typhoon. From August to October 1533 the Spanish Armada 
was beset by storm after storm which resulted in the sinking 
of many ships and failure of Spain's attempted invasion of 
Britain. On D-Day, June 6, 1 9 U- 4 , the allied invastion force 
took advantage of predicted good weather between storms and 
nade 4 -he amphibious assault across the beaches of Normandy. 
Tn December cf the same year Admiral Halsey's Third Fleet 
was caught preparing to refuel by an undetected typhoon, 
causing 28 ships to be crippled, 156 airplanes to lest, 
and the destroyers SFFNCS, MONAGHAN, and HULL to be sunk 
Wash 1576). 

Modern naval warfare technology is providing increas- 
inaly capable and complex, but environmentally dependent 
weapons/sensor systems. Net only the severe sea and weather 
conditions must be predicted to retain advantage over adver- 
saries, but also environmental factors which enhance or 
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degrade a wide range cf weapon and sensor systems which uti- 
lize electromagnetic (EM) propagation. Electromagnetic fre- 
quencies above the H? (3-30Mhz) band can be greatly affected 
by atmospheric refraction. 

a condition known as ducting occurs when a refractive 
layers cause EM energy to bend toward the earth at a rate 
greater than or equal to the earth’s curvature. Ducting 

occurs with certain critical vertical gradients of tempera- 
ture and humidity, and can cause both increased radar and 
radio ranges, and holes (gaps) in normal coverage. 

Tactical advantages exist by knowing duct locations. 
These advantages include being able to make realistic esti- 
mates cf ESM detection and counter detection ranges. This 
would enable commanders to make decisions concerning emis- 
sion control (EMCON) , the positioning of both air and sur- 
face surveillance assets, the altitude and flight profiles 
for strike aircraft to minimize detection, the placement cf 
-lectren jammers for maximum effect, and numerous other 
tactical considerations. 

Ducting commonly occurs with inversions which act as 
trapping layers, refracting or bending EM energy toward the 
earth. Inversions are stable layers between warm, dry air 
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Ducting is expected to be minimal near fronts and areas 
cf convective cloudiness. Fronts, with their associated 
upward motion, often dissipate the inversions as the whole 
air column becomes mixed. Areas of convective activity, 
astectible by the presence of cumulus clouds, are also nor- 
mally inversion free. Thus satellite infrared (IF) and vis- 
ible photography sho u.ld provide means to estimate ducting 
regimes. 

The Navy employs a microcomputer based system, IREPS 
(Integrated Refractive Effects Prediction System) (Hitnev, 
1979), “O Identify ducting conditions and to assess the 
effects or. various fleet EM emitters. IREPS requires radio- 
sonde data as input. 

The disadvantage of IREPS is that it predicts the duct- 
ing conditions only at the radiosonde launch site and only 
at the sounding time. As such, it is not an IREPS weakness 
but more an inherent weakness in single station assessments. 
The variability of the atmosphere makes projections in space 
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and predictions in tire of the i uctmg environment import ant 
for the tactician. while satellite data can provide some of 
the data necessary to make projections in space, something 



else is needed to predict chan 
EMCON, the ability tc forecast 
may be critical. Under strict 
cannot be launched because the 
mitting data could act like a b 

A predictive model exists 
atmospheric mixed-layer. The m 
developed by the Environmental 
Postgraduate School (Davidson e 
initialized with IREPS sounding 
of sea surface temperature, wi 
the inversion. 

Utiiizino the HP -884 5 (whit 
on board for use as part of th 
IREPS format, 3rower (198/) d? 
rate the NFS mixed layer mod? 
refract! vity profile for a 3 
reguires as inputs: surface wi 

and an estimated subsidence val 



es with time. Considering 
ducting for a 24 hour period 
EMCON conditions radiosondes 
ignal they emit while trans- 
aeon to hostile forces . 
for changes in the marine 
cro computer based model was 
Physics Croup at the Naval 
t al, 1933) . The model is 
data, surface observations 
nd speed, and subsidence at 

h ail carriers currently have 
e IREPS system) and data in 
v eloped a program to incorpo- 
1 to predict changes in ~ he 
3 hour period. The model 
nd, sea surface temperature, 
ue as well as current IREPS 
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sounding data to predict tie mixed layer evolution and 
resulting changes in both surface based and elevated duct 
heights up to 1200 m altitude. 

The purpose of this thesis is tc evaluate this state of 
the art atmospheric mixed layer model and ins application no 
fleet tactics. The data utilised in this study were from 
radiosondes taken by USS Niairz and USS C05AL SEA in the 
Arabian Sea in the vicinity of 2 3N 061E , a data-poor area, 
during February 1980. 
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II. THE MODEL DESCRIPTION 



A. MODEL OVERVIEW 

The NPS mixed layer model is an integrated, two layer, 
zero-order model which predicts changes i.i the marine atmos- 
pheric boundary layer (MABL) below the inversion. The two 
layers consist of the lower, well-mixed, turbulent boundary 
layer topped by the relatively non-i urb ul ent free a t m o ~ 
phere. The two layers are separated by ar. inversion or 
transition zone which is assumed to have zero thickness. 
The predicted properties are: 

(1) height of inversion 

(2) values of well-mixed properties 

(3) values of jumps at the inversion 

(4) formation of clouds/fog within the mixe 
Procedural uses are shown in Figure 2, and the egua 
will be presented in Section 3. Inputs to -he model are 

(1) Radiosonde data 

(a) Vertical distribution of temperature 

(b) Vertical distribution of moisture 

(2) Sea surface temperature 

(3) Wind speed 



:.a ye 



v./ u O 
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(4) Subsidence 

(5) Latitude 

(6) Julian date 

(7) Local time of sounding 

From the predicted properties, a refractivity profile is 
produced using *he refractivity equation (Section B) . 
Observed and predicted sea surface temperature and. vine! 
speed are input parameters. Subsidence rates are estimated 

ude, Julian date, and 
ider.t solar radiation. 



from previous observations. Latit 
local times are used to estimate inc 
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Figure 2 
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B. EASIC MODEL EQUATIONS 



Ones initial inputs have been made and the vertical 
structure digitized by fitting straight lines to approximate 
the vertical profiles of temperature and humidify, Equations 
1 A , IB, and 1C are used to calculate the flux scaling param- 
eters. The calculations are based on the air-sea tempera- 
ture and specific humidity differences and the wind speed. 
The scaling parameters are then used to estimate surface 
fluxes of momentum, heat, and moisture in Equations 2A, 23, 
and 2C. 

Equation 3A is used to predict the change in mixed-layer 
properties by the difference between surface and inversion 
fluxes scaled by inversion height. Equation 3B is used to 
predict changes in the inversion height by adding the flux 
at the inversion scaled by the inversion jump, a measure of 
the strength of the inversion, to the subsidence. Equation 
3C models the change in the inversion jump by combining the 



ef f 


ects of surface 


and inversion 


fluxes with 


subsidence and 




change cf the 


inversion he 


ight determine 


d in equation 


33. 


This requires 


use of the gr 


adient of the 


property above 


-he 


inversion. 
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Flux Scaling Parameters: 

u* = ca'^u 

I* = (9 - 9o) Ce' /2 

q* = Ce ' 2 (q - qc) 
Surface Fluxes: 

- (0 ' W ') = u * 2 

- (WTr* ) = □ *T * 

- ( rq ) = U *q * 

Where: 



(1A) 

(13) 

(1C) 

(2 A) 
(23) 
(2C) 



Cd,Ce = stability dependent drag coefficients 
U = wind speed 

= scaling parameters of vertical turbulent 
momentum, temp, and moisture transfer 
- (0 ' W ' ) = downward turbulent transfer of momentum 
- ( W ’ T ’ ) = upward turbulent transfer of heat 
- ( W 1 q 1 ) = upward turbulent transfer of moisture 



Predictive Equations (general form): 

dXm/dt = (WTx* c - W^X'i f f (?.) ) / h (3A) 

-WMC'i = DX (dh/dt - Ms) + f ( R) (33) 

d(DX)/dt = Yx ( dh/dt -Ws ) - (wHc' o-HMC' i+f (R) ) /h (3C) 
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Where: 



Xm = any well-ndxed property: 0, 9v, or g 

0 = potential temperature 
g = specific humidity 

W ' X * = turbulent vertical transfer (flux) of X 

0 subscript = "at the surface" 

1 subscript = "at the inversion" 

h = mixed-layer depth (or inversion height) 

DX = change in X across the inversion: (Xabove - Xm) 

Ws = subsidence at top of inversion 
Yx = vertical gradient of X above inversion 
f (R) = radiation factor 

Fefractivity Equation: 

M = 77.6 ?/T + 6 X 10 5 (q?/T 2 ) + . 1 57 2 (4) 

Where: 

M = modified r efr acrivit y 
? = atmospheric pressure in millibars 
T = temperature in Kelvins 
Z = altitude in meters 
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C. SATELLITE DATA INFUT TO MODEL 

Provision is made in the mod<=l for incorporating fore- 
cast wind and sea surface temperature changes during the 
prediction period. Although the mixed layer is sensitive to 
SST, in is believed that accuracies of */- 1°C in SST are 
sufficient for reasonable model accuracy. Stewart (1981) 
reports the following satellite measurement capability has 
been achieved: NOAA-6 (VHRR) infrared radiometer: SST to 

• 6°C (with r.c clouds) ; SeaSat and Nimbus-7 microwave radiom- 
eter (SUMP): SST to 1.0°C (with no rain, no RFI (radio fre- 

quency interference), and >600 km from land). 

Satellite SST could be used in the flux scaling equation 
(13) to determine 9c. Three methods are currently used 
aboard ship for SST determination: Bathythermographs, sea 

water injection temperature, and the ' bucket /thermometer ' 



method. None of 


the se 


yield the 


actual 


sea 


surface "skin” 


temperature which 


is the 


relevant 


quant it y 


2 . 


the aetermina- 


tion of stability 


and t 


urbuient 


heat ana 


mo 


isture fluxes. 



Therefore, satellite measurements would not only provide 
cloud descriptions over broad areas, but would provide a 
means to obtain the "skin" temperature. 
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III. ARABIAN SEA CLIMATOLOGY AND SYNOPTIC 
n|I SCRIPT TON 

A. CLIMATOLOGY 

The climatology of the northern Arabian Sea is dominated 
by the summer tsourhwest) and winter (northeast) monsoons. 
The data used in this study were from February during the 
latter parr of the northeast monsoon which is characterized 
by light offshore winds averaging 5-ID <nois, and normal 
Hadley circulation with the subtropical jet located near 30 
N at the 200 mb level. Fig. 4 shows the average per cent of 
time rhat elevated ducts occur during February to March in 
this region. 
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Figure 4, Frequency cf Elevated Duct. Occurrence Durinq 
Spring, after Lammers et al (1980). 



3. SYNOPTIC CONDITIONS 

FNCC (Fleet Numerical Oceanographic Center) products 
were used to describe the synoptic situation during the da^a 
period. They were: Surface Analyses, 500 mb analyses with 

Temperature fields, 300 mb analyses, 250 mb wind fields. 



Also, NOAA-6 


satellite 


photos w 


ere uss-i 


*■ O d 


e ter mine cloud 


coverage cf - 


he area as 


well as 


to ronfirn 


jet 


positions. 


The syr.cp 


tic situati 


o n in t h 


e region 


from 


6-23 February 



1980 was characterized by a stable high pressure surface 
ridge ever the Arabian peninsula and Arabian Sea, with low 
pressure troughs over the Red Sea to the west and Indian 
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subcontinent to the east. The most significant synoptic 
change occurred from 11 to 1 4 February as an 850 mb low 
moved across the northern part of the region. This low was 
northwest of the Persian Gulf a t 30 degrees North latitude 
on 11 February and tracked eastward, arriving over northern 
India on 14 February. It was during this brief period of 
11-14 February that the minimum shipboard surface pressures 
for the entire period were recorded. The minimum recorded 
surface pressure was 1012.0 mb on 14 February. 

Table I lists cloud coverage over rhe immediate area of 
the da~a . Low level stratus clouds were confirmed only on 
13 February and from 19 to 21 February. 



TA3LE I 

Data Area Cloud Coverage 9-21 February 1980. 



DATS CL COD TYPE PRESENT 



8-11 
1 1-12 

1 3 

14-18 

19-20 

2 0-21 

* aet ermin 



clear with scattered cirrus 

no satellite coverage 

low level scracus 

clear with sc at re red cirrus 

cirrus and stratus 

cirrus, stratus with multilevel 

ed from NO A A 5 IR imagery 



to 



NW 



An a t* emp* was made to derive subsidence values from 
synoptic scale wind fields by calculating divergence at var- 
ious levels. It was noted that at upper levels, wind maxima 
or jet screams had a dramatic influence or divergence pro- 
files. Although the magnitudes of vertical motion obtained 
in this manner were inconsistent and highly variable, an 
interesting end possibly significant observation was made. 

The subtropical jet associated with the wind maxima dis- 
cussed above is usually found at the northern limb of the 
Hadley Cell (Fig. 5). It is at this northern limb where the 
greatest downward motion or subsidence is found (Palmen and 
Newton, 1959). Although this is a feature of the large 

scale circulation, it seems evident '■hat the jet does have 
an effect which should be considered in single station 
assessments whenever possible. 

Changes in the subsidence rates calcalated by hindcast- 
ing and shown in Table II correlated closely with changes in 
the position and strength of the upper level jet scream. 
During the period 6-9 February there was a steady decrease 
in subsidence rates from -.0073 to -.0055 m/s. During the 
same period the tail of the 250 mb jet max moved steadily 
eastward from a point due north of the data area. From 6-11 
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February the 300 mb jet max shifted northward away from the 
data area. 



TAB L2 II 

Calculated Subsidence Fates 



RUN 


1 


2 3 


4 


5 5 


7 


8 


Q 


DATE 


6-7 


6-7 7-3 


8-9 


8-9 11-12 


14-15 


21-22 


22-23 


METHOD 


H 


-.30 


-.70 - . 65 


-.60 - 


.55 -.55 


-.75 


-.31 


-.50 


Q 


•f 


-1.72 + 


-.39 - 


.06 -1.17 


-.55 


-.41 


-. 33 


A 


-.54 


-.46 -1 . 10 

* subsidence 
** H hind cast 


+ 

rates 
, Q so? 


* -.37 

in o m /s. 
cific humi 


- .33 

ditv, A 


+ -.50 

adiabatic 



From 9-12 February, subsidence rates remained at -.0055 
m/s. Also, no change in strength or position of the 250 mb 
jet occurred her ween 11 and 12 February. 

A subsidence rare of -.0075 m/s was determined for the 
period between 14 and 15 February which was an increase 
since “he 12” h. The 300 mb jet strengthened and shifted to 
the south toward the data area from 12-15 February, and the 
250 mb jet increased in strength from 14-15 February. 

Curing the 21-23 February period, subsidence rates 
increased from -.0031 to -.3050 m/s. During this period. 
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bet h the 300 mb and the 250 mb jet oaxiaa moved closer :c 
the data area with the 300 mb is* strengthening and the 250 
mb jet weakening. 



Vertical 

Radiative convection 

sinking _ n ° 




Figure 5 . 



Hadley Cell and Subtropical Jet Position Shewing 
Area of Maximum Downward Mott on, after Palmer, and 
Newton ( 1 9 69) . 
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THE MODEL PERFORMANCE 



I V , 



A. DUCT PREDICTION 

Each sounding was run with I REPS to establish a baseline 
to assess the model's ability to predict changes in ducting 
conditions. It was -.her. possible -o ccmpare the model pre- 
dicted duct heights with those observe d. Seven different 
pairs of soundings were examined, with each pair separated 
by 24 hours. Two of the 24 hour periods had 2 simultaneous, 
spatially separated soundings at the end of the periods. On 
one occaisicn, one of the second soundings was taken at a 
location 425 nautical miles to the southeast of the first 
sounding . 

The observed winds at the start ana end of the period 
were used and interpolated linearly within the period to 
sp-cify the wind for the prediction period. Sea surface 
temperature at the latitude and longitude of the initial 
sounding and at the latitude and longitude of the verifying 
sounding were also interpolated within the period. 

Subsidence velocities were determined by hindcasting 
such that the predicted inversion height agreed with that 



28 



Predicted ducts using 



observed in the verifying sounding, 
rhis "hindcast " subsidence are compared with IREPS observed 
ducts in Table III . The SMS errors of the model are 
approximately one half the errors of persistence. Persis- 
tence is considered to be a "or edict ion" of no change. If 
is frequently used as a baseline against which predictions 
are compared to determine if the predictive method has 
merit. These results indicate that when subsidence can be 
accurately estimated, the model's performance is clearly 
superior to persistence in predicting due 1 -, bases and tops. 



TABLE III 

Comparison of Observed Duct Heiahts with Model Predicted 
Heights Using Hindcast Subsidence Values. 



RUN 


INITIAL 


VERIF YING 


MODEL 


NUMB 


ER SOUNDING 


SOUND ING 


PREDICTION 


1 


788-983 


4 70- 877 


345-740 


2 


788-983 


633- 938 


500-960 


3 


470- 877 


46- 504 


0-370 


4 


47-504 


0-135 


0-170 


5 


4 7-504 


0- 304 


0-156 


6 


no duct 


5-584 


1 1 5-375 


7 


~ 23- 1 03 1 
an d 0-112 


248- 569 


0-495 


8 


206-50 3 


no d uct 


50-420 


9 


no duct 


0- 258 


0-250 


RMS 


ERRORS: 


PREDICTION 


PERSISTENCE 




TOP 


173 


364 




BASE 


1 11 


214 




THICKNESS 


185 


285 



* heights in meters. 
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8 FEB OOZ 




Figure 6, Model Predicted Ducts (dashed) vs. Observed Duc~s 
(solid) for the Period 6-8 Fsoruary 1980. 

# 

For ? true predictive application, at least 3 sequential 
soundings are required. The first and second sounding 
enable the estimation of subsidence. The model is then ini- 
tiated from ".he second sounding using the subsidence 
obtained from the first and second. The third sounding is 

used to verify the prediction. Gleason (1982) found that a 
method which is based on the well mixed specific humidify is 
an accurate method to calculate subsidence. However, assum- 
ing persistence in subsidence rates appears to be the best. 
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a steady decrease in subsidence rates (Table II) during 
the first week of the data period coincided with a steady 
decrease in surface pressure. The upper level jet stream 
shifted away from the data area and weakened during this 
period. This suggests that judicious modification of oer- 
si stance -o btained subsidence may be appropriate in some 
cases especially when a trend has been observed. 

From three periods having 3 sequential soundings. Tables 
IV and V shew the results of predictions. In the first 
period (Run number 3) , an elevated layer became nearly sur- 
face Dasec . The model prediction was for a surface based 
duct. In the second period (Run number U) , a low elevated 

t 

duct became surface based. The model prediction matched the 
observation nearly perfectly. In the third case (Run number 
9) , a surface based duct formed when there was no duct ini- 
tially. In this case, the model predicted the formation of 
an elevated duct. While the predicted duct in the third 
case differed substantially in both heiaht and thickness 
from the observed duct, the model prediction was significant 
because the occurance of a duct was predicted from an ini- 
tial non -ducting condition. 



31 



For this particular comparison, the absolute errors ar 
more meaningful than rms errors. Two of Three model prs 
dice ions are essentially on the mark while all of the erro 
is concent ra red in the third. In the persistence cases, th 
errors are fairly evenly distributed. 



TABLE IV 

Comparison c” Observed Ducts with Model Predicted Ducts 
A ssuming Persiscent Subsidence. 



PUN 

NUMBER 



INITIAL 

SOUNDING 



VERIFYING 

SOUNDING 



MODEL 

PREDICTION 



3 

4 

Q 



470-877 
46-504 
no duct 



46- 504 
0 - 135 
0 - 258 



0-350 

0-140 

450-820 






he rg ire. s 



miners, 
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TABLE 7 



Error Analysis of Model. Prediction vs. Persistence. 



ABSOLUTE ERROR 



3 RUN RMS ERROR 



RUN 


MODEL 


PERSISTENCE 


MODEL 


PERSISTENCE 


3 TO D 


154 


373 


336 


337 


BASE 


46 


424 


284 


246 


THICKNESS 


108 


51 


75 


240 


4 TOP 


5 


369 






BASE 


0 


47 






THICKNESS 


5 


323 






9 TOP 


56 2 


253 






BASE 


49 0 


0 






THICKNESS 


72 


253 






* heights in met 


ers . 








CLOUD PREDICTION 










The cloud format i 


cn is 


predicted on the 


b as r s 


of cor.cur- 


prediction of th 


- lift 


ing condensation 


level 


(LCL) and 



the inversion height. If the LCL is below the inversion, 
clouds occur; if it ie above, no clouds occur. When model 
predictions showed no low level clouds forming, satellite 
imagery confirmed -.here were no low level clouds and than 
the sounding site either wets clear or showed scattered high 
cirrus clouds. The for mat i on of low level stratus clouds 
was predicted on 3 occasions over the entire period. Stra- 
tus clouds were predicted to form at 2100 local time 11 Feb- 
ruary. The first satellite imagery available following this 
time was at 0800 13 February, and showed the presence of low 
level stratus clouds (Fig. 8). 
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Figure 7. Model Prediction Sho 
Inversion at 2100 11 

Clouds were predicted to f 
hours into the 30 hour predict: 
dropped slightly below th e invs 
in the other cases. Satslli 
shows no clouds (Fig. 10) . 

Satellite imagery shows the 
tus clouds on one 2000 pas- 
amounts on the 0800 20 February 
were available to initialize *h 
formation of clouds, on the ne 
were predicted throughout the 
0700 22 February. Satellite im 



ing LC1 Droopir.a Beneath the 
FsB. 

rm at 0500 15 February (25 
or.). In this case the LCL 
sion, bur less sharply than 
e imagery from 15 February 



presence of 


low level srra- 


19 February 


and increased 


pass. Although no soundings 


model prior 


to this initial 


t available 


model run clouds 


period 0700 


21 February to 


gerv confirm 


the presence of 



3U 







Fiaure 8. NO A A 6 Satellite IS Imagery from 1800 13 February 
Showing the Presence of Stratus Clou is, 

both cirrus and stratus clouds on both the 0800 and 2000 
passes on 21 February (Fig. 12). 

Although the satellite imagery analysis was by neces- 
sity, subjective, some objectivity was retained by confirm- 
ing the initial analysis with a second analysis by an 
independent source. Both analyses agreed on cloud types. 
Since the model could be shown to be incorrect only on one 
occasion, the model’s overall performance in predicting the 
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MIXED LAYER (m) 




T!M£ 



Figure 9. Mcdsl Prediction 0303 1 4 "to 0900 15 February. 
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MIXED LAYER (rn) 




0700 1300 1 90C 0100 0700 1300 

TIME 



Figure 11. 



Model Prediction 0730 21 to 1300 22 February. 




Figure 12. NCAA 6 Satellite IP. Imagery 0300 21 February. 
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formation of low level clouds/ fog was sight our of r. in a 
cases correct or 893 correct., 

C. EXAMINATION CF MO EEL SENS IT I VIII 

Surface observation of sea surface temperature (SST) and 
wind can be obtained with an accuracy of +/- 1 °C and +/- 2 
knots re spec- ively . To test the model sensitivity to these 

possible measurement errors, successive model runs were made 
with SST 1 °C higher and lower chan observed. The effect of 
wind measurement error was evaluated using wind values 2 
knots higher ~har. observed. The resulting duct predictions 
were than compared with ducts predicted by actual observed 
winds and SSTs. Table VI contains the tabulated results. 

Noting that RMS error in duct specifications due to poten- 
tial SST measurement =rror are approximately double the RMS 
error due to potential wine, measurement error, one might 
reasonably conclude SST determination is the more critical 
of the two parameters. It is imperative to recognize that 
the wind has a much wider range of variability than the SST 
and accurate specification of both parameters are important 
for good model performance. 
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TABLE V! 



Sensitivity of Model Duct Predictions to SST and Wind 

V arratio ns. 



RUN 


VERIFYING 


SST 


SST 


WIND 


WIND 


NUMBER DUCTS 


+ 1°C 


-1°C 


+ 2 KTS 


-2 KTS 


1 


345- 740 


400-800 


3 15-7 15 


390-810 


330-725 


2 


500-960 


780- 1 190 


4 20-870 


625- 1030 


545-975 


3 


0-37 0 


0-395 


0-360 


0-400 


0-365 


4 


0-170 


0-270 


0-1 10 


0-180 


0- 140 


5 


0- 153 


0-210 


0-1 15 


0-170 


0-130 


6 


1 1 5-375 


170-570 


0-8 0 


25-310 


0-240 


7 


0-495 


290-430 


0-4 4 0 


0-5 10 


0-450 


8 


50-420 


0-360 


0-3 10 


70-450 


0-370 


9 


0- 25 0 


0-240 


0-290 


0-240 


0-190 


Combined 


RMS error: 


SST 


DU CT 


» I N D 








114 


top 


5 5 








104 


base 


5 0 








107 


hi ckne ss 


2 6 




* ha 


ights in mete 


rs. 








The 


model's demonstrated 


s er.si ti v 


ity to sea surface 


temperature points to 


^ h - ri 0 s 


1 for accurate SST 


det ermina- 


tions. 


Current methods to m 


s a sure sa 


a surface t 


emperature 


in cl ude: 


mcr.it or ir. g 


sea ware 


r injecti 


on temperat 


ure (20-40 


ft. ben 


eath the surf 


ace) ; ut 


ilizinc E7 drops (ge 


t tempera- 


fure at 


afcou~ 1 ft . ) ; 


and the 


bucket a 


nd thermometer method 


(get temperature of top few in 


oh es) 


None of the 


ss aietnc as 



yield the sea surface skin n9.Tiparat.ure which is the quantify 
necessary fc determine the surface fluxes. These fluxes, in 
turn, drive the evolution of the mixed layer, the inversion, 
and thus ducting conditions of the lower atmosphere. 
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The 052340Z February 1980 sounding was used to test 
model sensitivity to varying subsidence rates. The hindcast 
subsidence was -.0030 m/s (Table II) which resulted in a 
predicted 355 meter thick elevated dues with the top an 740 
meters elevation and the base at 345 meters. Subsidence was 
varied from this baseline value and deviations in the pre- 
dicted duct height and thickness were recorded. This par- 
ticular sounding was chosen because the associated initial 
and verifying ducts were elevated, which allowed the varied 
subsidence tc affect the position of the predicted duct in 
both the upward and downward directions. 

Deviations in the duct predictions resulting from subsi- 
dence specification errors of +/- .0030 m/s are 245 m at the 
top, 269 m at the base, and 25 m in thickness (Table VII). 
Errors resulting from assuming persistence in duct heights 
(Table III) are approximately equal to the errors resulting 
from a subsidence specification error of -/- .0030 m/s. 
From this single case comparison, it is clear that the mo del 
is an improvement 2Y.iT Eiisistence only when subsidence can 
be specified to wpthin ♦ /- • 003 0 m/s . 

The range of subsidence calculated during the short 
period from 6 to 22 February 1980 in the data area was 
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-.0030 to -.0075 m/s or -.00525 +/- .00225 m/s. Assuming 

this to be a normal range of variability and the median 
value, -.00525 m/s, tc be representative of the climatologi- 
cal subsidence, then the "climatological’" value appears to 
be less than the +/-.0030 m/s marginal utility criterion. 
This conclusion, based or. a limited data set, seems to indi- 
cate that a climatological subsidence value would provide a 
duct height forecast which, on the average, has some shill 
over persistence or a forecast of no change. 

In reali-y, errors will not be confined to subsidence 
alone, but will be caused by errors in other inputs such as 



wind and SST. While 


SST, wind. 


and subs id 


er. ce 


errors may be 


offsetting at times. 


there will 


be times 


when 


they will he 


additive. When all 


the errors 


are addit 


i ve , 


the expected 



duct height error for SST plus wind measurement errors are 
169 m at the top and 154 tn at the bottom {Table V) . Given 
the errors due to assuming persistence in duct heights of 
364 m at the top ana 214 m at the base of the duct (Table 
IV) , the difference of approximately 203 m at the ~cp and 60 
m at the base indicate the maximum duct height error permis- 
sible due to subsidence specification error is, in the mean, 
130 m. From Table VI, the subsidence specification error 
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expected to produce this 130 m duct height error is approxi- 
mately + /- .0015 m/s. 

From the cases (Table III) in which the best possible 
estimate of subsidence, the hindcast value, was used, the 
FMS model prediction error averages 143 m less than persis- 
tence error. I~ will be shown later that if subsidence is 
allowed to vary from this "best estimate", a deviation of 
about +/- .0017 m/s will produce the additional 143 m error 
in height of duct top and base necessary to render the pre- 
diction no better than persistence. 

In conclusion, evaluations of the model with the data 
indicate that with the expected SST and wind measurement 
errors, subsidence must be specified to within +/- .0015 m/s 
for the model to consistently outperform persistence. The 
observed range of subsidence variablity from this limited 
data set of +/- .00225 m/s effectively rules out the use of 
a climatological subsidence as a viable input parameter. 
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TABLE VII 

Sensitivity of Model Duct Predictions to Subsidence. 



SUBSIDENCE 




PREDICTS 


D DUCT PARAMETERS 




(m/s) 


TOP 


CH ANGE 


BASE 


CHANGE 


THICKNESS 


CHANG 


.0000 1 


020 


2 80 




645 


300 


375 


-20 


-.0005 


970 


2 30 




585 


240 


385 


-10 


-.0010 


9 10 


1 70 




535 


1 90 


375 


-20 


-.0015 


865 


1 25 




480 


135 


385 


-10 


-.00 20 


830 


90 




430 


85 


400 


5 


-.00 25 


7 8 5 


45 




380 


35 


405 


10 


-.0030 * 


740 


0 




34 5 


0 


395 


0 


-.0035 


7 1 3 


-27 




312 


-33 


401 


6 


-.0040 


675 


- 65 




26 0 


-85 


415 


20 


-.0045 


630 


-1 10 




225 


-120 


405 


10 


-.0050 


600 


- 1 40 




180 


-165 


420 


25 


-.0055 


570 


-1 70 




145 


-200 


425 


30 


-.0060 


535 


-2 05 




110 


-235 


425 


30 


EMS ZEROES 


z 














DEVI 


ATION FFOM 


REFERENCE 


EXPECTED DEVIATION 


SUBSIDENCE 


TOP BASE 


TH ICKNESS 


(Ws X 24 


hours) 


+/- .0030 


245 


269 




25 


259 




+/- .0025 




202 


221 




22 


216 




+/- .0020 




1 56 


178 




23 


173 




+/- .0015 




1 18 


1 28 




10 


130 




+/- .0010 




79 


85 




15 


86 




+/- .0005 




37 


34 




8 


43 





* reference values 
** heights in meters 
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V. 



TACTICAL APPLICATIONS 



The first consideration in tactical application of the 
model is the recognition of its capabilities and limita- 
tions. It can be used to estimate when low clouds or fog 

can be expected to form. It can be used to estimate the 
expected position of low elevated and surface based ducts 
below 1200 m over an 18 - 24 hour period. However, it can- 
not predict information on upper level ducts. 

The tactical significance of knowing up to 30 hours in 
advance that visibility may be reduced or a low ceiling 
developed will depend on current and planned operations. 

t 

The results presented in Chapter IV B indicate that the 
model does well in predicting low level stratus/fog. 

When fcq or low level stratus clouds are predicted, the 
following possible effects on own and enemy forces in each 
warfare area should be considered: 

(1) maneuverability restricted 

(2) lest or reduced visual signalling capability 

(3) reduced visual target detection/identification 

(4) degraded optical equipment performance 
including 13 weapons and sensors 

(5) hampered flight operations for carrier based 
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aircraft and LAMPS 



Tactical employment of Naval assets requires a knowledge 
of the state of the environment if those assets are to be 
employed effectively. To that end, IREPS (Integrated 
Refractive Effects Prediction System) was developed. when 
radiosonde data are inpur. to IREPS, a r ef ractivi^.y profile 
is generated which is utilized to define the location of 
ducts and to assess radar coverage, ESM and communication 
ranges. 

Because soundings are normally taken only once or twice 
daily from aircraft carriers and the atmospheric boundary 
.iyei undergoes change, rhere exists a need to be able *o 
predicr changes in the refract! vity profile between sound- 
ings. This prediction can be critically important to battle 
croups in the positioning of both surface and air assets, 
ana in the intelligent management of the 3MC0N plan. 

The duct prediction capability, as detailed in Chapter 
Z'J A, provides a significant planning tool. Tactical useage 
of the duct predictions stems from the effects of ducts on 
EM propagation. Prediction of a surface based duct indi- 
cates probable extended ranaes for transmitter-receiver 
antenna pairs in the duct and possible holes or gaps in 
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near, or below 



coverage just above the duct for antennas in, 
the duct . 

The predicted M profile may be used as an environmental 
data set. to run IREPS and coverage diagrams generated for 
specific emitrers. This procedure was used with a hypothet- 
ical surface search radar under observed refractive condi- 
tions and model predicted refractive conditions. IREPS 
version 1*7 (unclassified) was used to generate the coverage 
diagram 1 ;. Radar parameters for the hypothetical radar were: 
Antenna height 110 ft., antenna type sin(x)/x, vertical beam 
width 10 ° r elevation angle 0°, free space range 55 nautical 
miles, and frequency 5 Ghz. 

Cleon.'. y shown in Figs. 14 and 15 is the model's value in 
predic-f r.q the occurance of low level aucts. The radar 
coverage an Fig. 14 is significant in that the initial pro- 
pagation conditions were for surface target detections out 
to about 25 nin and the predicted conditions were for detec- 
tions well in excess of 100 nm. The predicted propagation 
conditions were in close agreement with those observed. 

Ar. example of how the model might be used operationally 
is as follows: 

(1) radiosonde is launched and surface observations 
of SST and wind are recorded. 
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HEIGHT (ft) 





Figure 14. Surface Search Radar Coverage 0333 8 Feb 80 
Based on Actual Sounding (1) and Model 
Prediction (r) . 
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HEIGHT (ft) 




Figure 15. Surface Search Radar Coverage 3240 9 Feb 80 
Based on Actual Sounding (1| and Model 
Prediction (r) .. 

(2) radiosonde data are input to IREPS and coverage 
diagrams are Generated to make current assessments. 

(3) radiosonde data and surface observations are input 
to the NPS model and a 30 hour prediction is mads. 



(4) -.he predicted 1 profile is input to IREPS and 
predicted coverage diagrams are generated. 

(5) operational plans, incorporating forecast ducting 
effects, are proposed. 

(6) tactical decisions, using the environment to the 
advantage of the force, are made. 
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In a scenario similar to that observed in the Arabian 
Sea 7 to 9 February 1980 when an elevated duct became sur- 
face based, many tactical considerations are involved. Sup- 
pose the bar. tie group wanted to avoid detection by an enemy 
force in the region. Forewarned that surface based ducting 
was expected to occur in 12 to 18 hours and that OHF radio 
transmissions,, normally useful between ships only up to 25 
to 30 nautical miles, could be defected at ranges in excess 
of 200 nautical miles due to ducting, the battle group com- 
mander would be able to impose a more restrictive EMCON con- 
dition in adequate time to prevent hostile intercept of the 
group’s transmissions. Alternatively, if the battle group 
commander expected to be opposed by active surface search 
radar, he could: initiate a preemptive attack before his 
presence was exposed; steam out of range; or, relax EMC ON at 
the appropriate time. 

For a battle group in a multithreat environment, early 
threat detection is critical to defeat the threat at minimum 
cost to the battle group's offensive capability. The 12 to 
18 hour advance notice of surface based ducting provides the 
tactical planner the -ool necessary tc use his limited sur- 
veillance assets with the efficiency needed for sustained 
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operations and the e ffectivenes s require! for early threat 
detection. Surface ships could be positioned in such a way 
as to take maximum advantage of the extended surface search 
radar ranges afforded by the surface based duct. Airborne 

surveillance assets,, freed from task of surface surveil- 
lance, could concentrate their efforts toward detecting and 
tracking air targets The net result is increased force 
ef f ectiver.ess . Depending on the battle group's speed and 

the picket station assigned, th® surface ship may require 
several hours to reach the assigned station, thus making the 
lead time provided by the prediction all the more important. 

The timing of war at sea strikes to coincide with pre- 
dicted occurrence of t surface based duct could give consid- 
erable afivan-.aae to the strike aircraft by enabling them to 
use predicted holes jus t above the duct to penetrate enemy 
surveillance and to us > the duct for targeting and jamming. 
In such a case, a 12 to 13 hour prediction provides time to 
plan mission tactics * arm aircraft, and conduct mission 
briefs. 

When sonobuoy patterns are laid and monitored by LAMPS, 
one constraint on placement of the bouys is LAMPS' on sta- 
tion endurance. When a surface based duct occurs, sonobouys 
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can be monitored at greater distances, 
advance of surface based duct in g provides 
essary to plan and lay bouv patterns at 
battle group which take advantage of the s 
range. 



Prediction in 
the lead time nec- 
di stances from the 
xtended monitoring 
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VI 



• CONCLUSIONS AND REC0HJJ3NDA7I0NS 

The Tactical Environmental Support System (TESS) , cur- 
rently under development, will incorporate various satellite 
derived products, and single station assessment systems such 
as IREPS and the NPS mixed layer model. This will place 
state of the art environmental sensing and predictive capa- 
bility at the disposal of those who need the information: 
the operational commanders. 

Eased on evaluations with a limited data set, the NPS 
model has been demonstrated -o perform well both in pre- 
dictions of lew level cloud formation and in predictions of 
low elevated and surface based ducts. Further work is 
required to determine the modal’s range of applicability 
both geographically and seasonally. It Ls believed that 
when subsidence can be specified with an accuracy of 
•*•/-. 0015 m/s, the model can be used to accurately forecast 
low level ducts. Two additional recommendations are 
detailed in the following paragraphs. 

In tactical situations where the outcome may depend on 
accuracy of duct predictions, a method is needed to verify 
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the prediction. A radiosonde launch may be impractical due 
to EMCON conditions. An acoustic sounder can be used to 
measure the height of the inversion which, as shown in Fig. 
3, corresponds -o the top of the duct. Davidson , et al 
(1982) , while conducting marine atmospheric boundary layer 
research, utilized the economical acoustic sounder on R. V. 
ACANIA to monitor inversion height. Using a combination of 
daily radiosondes, NPS model duct predictions, and periodic 
acoustic soundings to track the inversion height would pro- 
vide improved caoability to use the environment. It is 
therefore recommended that similar acoustic sounders be 
installed on aircraft carriers. 

Determining appropriate temperature ana humidity pro- 
files during initialization of the monel is a possible 
source of error. There are two schools of thought which 
differ on the appropriate height fcr the inversion initiali- 
zation. Sounding data typically show the inversion as a 
layer 100 to 200 meters thick, while the model requires a 
single height to be specified. One school suggests that the 
bottom of the inversion layer is best while ■‘■he second sug- 
gests that the midpoint of the layer is most reasonable. 
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The work in this thesis used the forier method. The 
model predictions (Table III) show a bias in that nearly all 
of the predicted duct heights are lower than those observed. 



This seems 


to indicate that 


i nver sion 


heights should 


b 9 


digitized at 


the midpoint 


of 


the layer 


instead of at 


t he 


bottom. It 


is recommended 


that 


further w 


ork be done in t 


his 



area to establish a mere objective method of digitization. 

The NPS atmospheric boundary layer model could provide 
the fleet a predictive capability and '■acticai planning 
fool. Consideration should be given to introducing it into 
the fleet for operational evaluation. The need to under- 

stand the model assumptions in order to assess model apoii- t 
cability in a given environmental situation, the need for 
accurate subsidence estimates, and the requirement for sur- 
face wind predictions preclude model installation except 
where qualified personnel are assigned. In this regard, the 
model should be installed at locations (aircraft carriers) 
where there are Geophysics Officers who can interpret and 
use the model as a tool along with other meteorological 
information. 
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APPENDIX A 



RADIOSONDE TINES AND LOCATIONS 



TABLE VIII 
Radiosonde Data Set. 



SOUNDING 


SHIP 


DTG 


LATITUD 


1 


NIMITZ 


0523 40Z 


2 3 0 1 2 1 N 


2 


MIMITZ 


0700 20Z 


22 ° 3 8 * N 


3 


CORAL SEA 


070030Z 


1 7 °40 ' N 


4 


NI MITZ 


0723 33Z 


23° 0 7* N 


5 


CORAL SEA 


082240Z 


23°00 *N 


6 


NI MITZ 


0823 33Z 


23 0 1 3 • N 


7 


NIMITZ 


1 100 1 4 Z 


23° 1 7* N 


8 


NIMITZ 


1 123 46 Z 


23 o 1 0 ' N 


9 


NI MITZ 


1 323 07Z 


23°03 *N 


10 


NIMITZ 


1500 1 3 Z 


23° 1 1 » N 


1 1 


NI MITZ 


2103 00Z 


24 ° 0 1 ' N 


12 


CORAL SEA 


2120 35Z 


23 ° 0 6 ’ N 


13 


NIMITZ 


212244Z 


23°13' N 


14 

■JU 


CORAL SEA 

r ■ _ _ .. . 


2220 1 5 Z 


23 ° 3 8' N 



* date-time groups February 1930 



TABLE IX 

Radiosondes Used in Model Performance Runs 



RUN NUMBER 1 

SOUNDINGS : 

INITIAL 1 

VERIFYING 2 



2 3 4 5 

12 4 4 

3 4 5 6 



6 7 

7 9 

8 10 



LONGITUD 



061° 09 
061^36 
06 5 °5 1 
061000 
061°11 
061° 00 
060° 5 9 
06 1 °0 1 
064° 30 
06 1 °07 

05 8 ° 3 2 
061°12 

06 0 °5 8 
060 o 0 9 



8 

1 1 
12 



1.TJ hr.l t*j til L X J tiJ W w tu W DJ til td 
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